home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 15 / CU Amiga Magazine's Super CD-ROM 15 (1997)(EMAP Images)(GB)[!][issue 1997-10].iso / CUCD / Graphics / Ghostscript / source / zarith.c < prev    next >
C/C++ Source or Header  |  1997-03-04  |  7KB  |  292 lines
  1. /* Copyright (C) 1989, 1992, 1993, 1994, 1997 Aladdin Enterprises.  All rights reserved.
  2.   
  3.   This file is part of Aladdin Ghostscript.
  4.   
  5.   Aladdin Ghostscript is distributed with NO WARRANTY OF ANY KIND.  No author
  6.   or distributor accepts any responsibility for the consequences of using it,
  7.   or for whether it serves any particular purpose or works at all, unless he
  8.   or she says so in writing.  Refer to the Aladdin Ghostscript Free Public
  9.   License (the "License") for full details.
  10.   
  11.   Every copy of Aladdin Ghostscript must include a copy of the License,
  12.   normally in a plain ASCII text file named PUBLIC.  The License grants you
  13.   the right to copy, modify and redistribute Aladdin Ghostscript, but only
  14.   under certain conditions described in the License.  Among other things, the
  15.   License requires that the copyright notice and this notice be preserved on
  16.   all copies.
  17. */
  18.  
  19. /* zarith.c */
  20. /* Arithmetic operators */
  21. #include "math_.h"
  22. #include "ghost.h"
  23. #include "errors.h"
  24. #include "oper.h"
  25. #include "store.h"
  26.  
  27. /****** NOTE: none of the arithmetic operators  ******/
  28. /****** currently check for floating exceptions ******/
  29.  
  30. /* Define max and min values for what will fit in value.intval. */
  31. #define min_intval min_long
  32. #define max_intval max_long
  33. #define max_half_intval ((1 << (size_of(long) / 2 - 1)) - 1)
  34.  
  35. /* Macros for generating non-integer cases for arithmetic operations. */
  36. /* 'frob' is one of the arithmetic operators, +, -, or *. */
  37. #define non_int_cases(frob,frob_equals)\
  38.  switch ( r_type(op) ) {\
  39.   default: return_op_typecheck(op);\
  40.   case t_real: switch ( r_type(op - 1) ) {\
  41.    default: return_op_typecheck(op - 1);\
  42.    case t_real: op[-1].value.realval frob_equals op->value.realval; break;\
  43.    case t_integer: make_real(op - 1, (double)op[-1].value.intval frob op->value.realval);\
  44.   } break;\
  45.   case t_integer: switch ( r_type(op - 1) ) {\
  46.    default: return_op_typecheck(op - 1);\
  47.    case t_real: op[-1].value.realval frob_equals (double)op->value.intval; break;\
  48.    case t_integer:
  49. #define end_cases()\
  50.   } }
  51.  
  52. /* <num1> <num2> add <sum> */
  53. /* We make this into a separate procedure because */
  54. /* the interpreter will almost always call it directly. */
  55. int
  56. zop_add(register os_ptr op)
  57. {    non_int_cases(+, +=)
  58.        {    long int2 = op->value.intval;
  59.         if ( ((op[-1].value.intval += int2) ^ int2) < 0 &&
  60.              ((op[-1].value.intval - int2) ^ int2) >= 0
  61.            )
  62.            {    /* Overflow, convert to real */
  63.             make_real(op - 1, (float)(op[-1].value.intval - int2) + int2);
  64.            }
  65.        }
  66.     end_cases()
  67.     return 0;
  68. }
  69. int
  70. zadd(os_ptr op)
  71. {    int code = zop_add(op);
  72.     if ( code == 0 ) { pop(1); }
  73.     return code;
  74. }
  75.  
  76. /* <num1> <num2> div <real_quotient> */
  77. private int
  78. zdiv(register os_ptr op)
  79. {    register os_ptr op1 = op - 1;
  80.     /* We can't use the non_int_cases macro, */
  81.     /* because we have to check explicitly for op == 0. */
  82.     switch ( r_type(op) )
  83.        {
  84.     default:
  85.         return_op_typecheck(op);
  86.     case t_real:
  87.         if ( op->value.realval == 0 )
  88.             return_error(e_undefinedresult);
  89.         switch ( r_type(op1) )
  90.            {
  91.         default:
  92.             return_op_typecheck(op1);
  93.         case t_real:
  94.             op1->value.realval /= op->value.realval;
  95.             break;
  96.         case t_integer:
  97.             make_real(op1, (double)op1->value.intval / op->value.realval);
  98.            }
  99.         break;
  100.     case t_integer:
  101.         if ( op->value.intval == 0 )
  102.             return_error(e_undefinedresult);
  103.         switch ( r_type(op1) )
  104.            {
  105.         default:
  106.             return_op_typecheck(op1);
  107.         case t_real:
  108.             op1->value.realval /= (double)op->value.intval; break;
  109.         case t_integer:
  110.             make_real(op1, (double)op1->value.intval / (double)op->value.intval);
  111.            }
  112.        }
  113.     pop(1);
  114.     return 0;
  115. }
  116.  
  117. /* <num1> <num2> mul <product> */
  118. private int
  119. zmul(register os_ptr op)
  120. {    non_int_cases(*, *=)
  121.        {    long int1 = op[-1].value.intval;
  122.         long int2 = op->value.intval;
  123.         long abs1 = (int1 >= 0 ? int1 : - int1);
  124.         long abs2 = (int2 >= 0 ? int2 : - int2);
  125.         float fprod;
  126.         if (    (abs1 > max_half_intval || abs2 > max_half_intval) &&
  127.             /* At least one of the operands is very large. */
  128.             /* Check for integer overflow. */
  129.             abs1 != 0 &&
  130.             abs2 > max_intval / abs1 &&
  131.             /* Check for the boundary case */
  132.             (fprod = (float)int1 * int2,
  133.              (int1 * int2 != min_intval ||
  134.              fprod != (float)min_intval))
  135.            )
  136.             make_real(op - 1, fprod);
  137.         else
  138.             op[-1].value.intval = int1 * int2;
  139.        }
  140.     end_cases()
  141.     pop(1);
  142.     return 0;
  143. }
  144.  
  145. /* <num1> <num2> sub <difference> */
  146. /* We make this into a separate procedure because */
  147. /* the interpreter will almost always call it directly. */
  148. int
  149. zop_sub(register os_ptr op)
  150. {    non_int_cases(-, -=)
  151.        {    long int1 = op[-1].value.intval;
  152.         if ( (int1 ^ (op[-1].value.intval = int1 - op->value.intval)) < 0 &&
  153.              (int1 ^ op->value.intval) < 0
  154.            )
  155.            {    /* Overflow, convert to real */
  156.             make_real(op - 1, (float)int1 - op->value.intval);
  157.            }
  158.        }
  159.     end_cases()
  160.     return 0;
  161. }
  162. int
  163. zsub(os_ptr op)
  164. {    int code = zop_sub(op);
  165.     if ( code == 0 ) { pop(1); }
  166.     return code;
  167. }
  168.  
  169. /* <num1> <num2> idiv <int_quotient> */
  170. private int
  171. zidiv(register os_ptr op)
  172. {    register os_ptr op1 = op - 1;
  173.     check_type(*op, t_integer);
  174.     check_type(*op1, t_integer);
  175.     if ( op->value.intval == 0 )
  176.       return_error(e_undefinedresult);
  177.     if ( (op1->value.intval /= op->value.intval) ==
  178.         min_intval && op->value.intval == -1
  179.        )
  180.        {    /* Anomalous boundary case, fail. */
  181.         return_error(e_rangecheck);
  182.        }
  183.     pop(1);
  184.     return 0;
  185. }
  186.  
  187. /* <int1> <int2> mod <remainder> */
  188. private int
  189. zmod(register os_ptr op)
  190. {    check_type(*op, t_integer);
  191.     check_type(op[-1], t_integer);
  192.     if ( op->value.intval == 0 )
  193.         return_error(e_undefinedresult);
  194.     op[-1].value.intval %= op->value.intval;
  195.     pop(1);
  196.     return 0;
  197. }
  198.  
  199. /* <num1> neg <num2> */
  200. private int
  201. zneg(register os_ptr op)
  202. {    switch ( r_type(op) )
  203.        {
  204.     default:
  205.         return_op_typecheck(op);
  206.     case t_real:
  207.         op->value.realval = -op->value.realval;
  208.         break;
  209.     case t_integer:
  210.         if ( op->value.intval == min_intval )
  211.             make_real(op, -(float)min_intval);
  212.         else
  213.             op->value.intval = -op->value.intval;
  214.        }
  215.     return 0;
  216. }
  217.  
  218. /* <num1> ceiling <num2> */
  219. private int
  220. zceiling(register os_ptr op)
  221. {    switch ( r_type(op) )
  222.        {
  223.     default:
  224.         return_op_typecheck(op);
  225.     case t_real:
  226.         op->value.realval = ceil(op->value.realval);
  227.     case t_integer: ;
  228.        }
  229.     return 0;
  230. }
  231.  
  232. /* <num1> floor <num2> */
  233. private int
  234. zfloor(register os_ptr op)
  235. {    switch ( r_type(op) )
  236.        {
  237.     default:
  238.         return_op_typecheck(op);
  239.     case t_real:
  240.         op->value.realval = floor(op->value.realval);
  241.     case t_integer: ;
  242.        }
  243.     return 0;
  244. }
  245.  
  246. /* <num1> round <num2> */
  247. private int
  248. zround(register os_ptr op)
  249. {    switch ( r_type(op) )
  250.        {
  251.     default:
  252.         return_op_typecheck(op);
  253.     case t_real:
  254.         op->value.realval = floor(op->value.realval + 0.5);
  255.     case t_integer: ;
  256.        }
  257.     return 0;
  258. }
  259.  
  260. /* <num1> truncate <num2> */
  261. private int
  262. ztruncate(register os_ptr op)
  263. {    switch ( r_type(op) )
  264.        {
  265.     default:
  266.         return_op_typecheck(op);
  267.     case t_real:
  268.         op->value.realval =
  269.             (op->value.realval < 0.0 ?
  270.                 ceil(op->value.realval) :
  271.                 floor(op->value.realval));
  272.     case t_integer: ;
  273.        }
  274.     return 0;
  275. }
  276.  
  277. /* ------ Initialization table ------ */
  278.  
  279. BEGIN_OP_DEFS(zarith_op_defs) {
  280.     {"2add", zadd},
  281.     {"1ceiling", zceiling},
  282.     {"2div", zdiv},
  283.     {"2idiv", zidiv},
  284.     {"1floor", zfloor},
  285.     {"2mod", zmod},
  286.     {"2mul", zmul},
  287.     {"1neg", zneg},
  288.     {"1round", zround},
  289.     {"2sub", zsub},
  290.     {"1truncate", ztruncate},
  291. END_OP_DEFS(0) }
  292.